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Abstract: Nuclear modification factors of inclusive prompt photon production in d-Au 
collisions at RHIC and p-Ph collisions at the LHC are provided at different rapidities. 
The calculations are performed at NLO accuracy using the EPS09 NLO nuclear parton 
distribution functions (nPDFs) and their error sets. The results are compared to the ones 
obtained with the nDS and HKN07 NLO nPDFs, and to the corresponding nuclear modifi- 
cation factors of neutral pion production in these collisions. The sensitivity of these results 
to the scale choice is also investigated. Interestingly, the predictions using the different 
nPDF sets differ from each other to the extent that this observable can be expected to 
become very useful for probing nPDFs over a wide range of Bjorken-x. In order to obtain a 
perturbative QCD baseline in heavy-ion collisions, calculations are carried out for minimum 
bias Au-Au collisions at RHIC and Pb-Pb collisions at the LHC. We also estimate the 
maximal possible suppression which the produced QCD matter can be expected to have on 
inclusive prompt photon production due to the quenching of the fragmentation component. 
The nuclear modification factor for prompt photon production is thus suggested to be used 
for gauging both the cold and the hot nuclear matter effects on other hard processes which 
are expected to be affected by quark-gluon plasma formation, such as large-p^ hadron and 
jet production. 
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1 Introduction 

The production of inclusive and isolated prompt photons in hadronic collisions is a process 
that is known - along with jet production - to carry direct information on the gluon Parton 
Distributions Functions (PDF) in the proton [1, 2]. However, prompt photon production 
data are presently not used in the global fits of PDFs but rather provide an independent 
test of perturbative QCD and of the universality of the partonic densities. The good 
overall agreement between next-to-leading order (NLO) QCD calculations [3] and world 
data from fixed-target (y/s = 23 GeV) to collider experiments {y/s = 1.8 TeV) indicates 
that the underlying parton dynamics is properly understood. In this respect, the expected 
high-precision measurements to be performed at the LHC and RHIC might give extra 
constraints on the PDFs in the gluon sector [4]. The agreement between the first ATLAS 
and CMS measurements [5, 6] and NLO calculations is thus one more encouraging step 
towards this goal. 

The recent global analyses of nuclear PDFs (nPDF) [7-14] have shown that nuclear 
collisions seem to follow the same collinear factorization theorem that works in hadronic 
collisions: Structure function data in nuclear (lepton) deep inelastic scattering (DIS) as 
well as Drell-Yan dilepton (DY) and hadron production measurements in high-energy 
proton(deuterium)-nucleus {p{d)-A) collisions are consistent with universal (process in- 
dependent) nuclear modifications of PDFs [12, 13]. Other checks have been carried out 
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recently using neutrino-DIS data [15, 16]. Whether such a good description persists for 
prompt photon production has however not been verified yet, as no decisive-precision p-A 
or d-A data are presently available. Keeping in mind the ongoing efforts at RHIC and 
appreciating the potential capabilities of the LHC, this situation is expected to undergo a 
change in the near future. 

As discussed in [17, 18], prompt photon production in p-A collisions appears to be a 
promising tool for probing the gluon nuclear densities, which up to now have been con- 
strained through scaling violations of the nuclear DIS and DY data and inclusive vr'^ pro- 
duction in d-Au collisions at RHIC [13] but, clearly, for which more constraints would 
be badly needed. Therefore, complementary to recent studies of inclusive pion produc- 
tion [19] and prompt photon+heavy quark production [20], we present and discuss in the 
present paper the NLO predictions for the nuclear modification ratios of prompt photon 
spectra in d-Au collisions at RHIC and p-Ph collisions at the LHC using modern nPDF 
sets (section 3). Baseline predictions in heavy-ion {A- A) collisions will also be discussed 
in detail in section 4. Before this, let us first briefly look at the theoretical framework used 
in the present analysis. 

2 Theoretical framework 
2.1 Elements of the calculation 

We consider in this analysis inclusive production of prompt photons, 

/ii + /i2 ^ 7 + X, 

in d-Au and Au-Au collisions at RHIC {^/s^ = 200 GeV) as well as in p-Ph and Pb- 
Pb collisions at the LHC (y^s^ = 8.8 TeV and ^s^^ = 5.5 TeV, respectively) nominal 
energies. 

Using collinear factorization, the inclusive photon production cross section can be 
written convolution 

dah,h,^y+x = Y,f-'{M^)^fj'{M^)^daij^^+x'{fi^M^) + 0{l/M^), (2.1) 

where (8) is an integral over the longitudinal momentum fractions xi and X2 of the incoming 
partons i and j, whose number densities - the PDFs of the projectile and target PDFs 
- are f!^^{xi, M^) and fj^{x2,M'^). The factorization scale entering the PDFs and the 
perturbatively calculable partonic pieces (subcross-sections at LO) daij^^+x' is denoted 
by M, and the renormalization scale by As usual, X indicates the inclusive nature of 
the cross section and X' indicates that in each partonic hard process (2-to-2 or 2-to-3) we 
integrate over everything else but the photon. 
At leading order (LO) in pQCD, 0{a 

cmC^s)) photons can be produced directly via 
two types of partonic subprocesses: Compton scattering q{q)g — t- "yq{q) and annihilation 
qq — )• jg. In terms of the transverse momentum p^, and rapidity y of the photon, and 
the cms-energy ^/s, the momentum fractions typically probed by direct photon production 
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are xi^2 ~ (2Pr/\/s)e^^. Thus, towards smaller p^, larger y and larger ^/s, the process 
becomes sensitive to PDFs at smaller X2 (target) and larger xi (projectile). At small 
values of x, the gluon distribution is larger than that of sea quarks, which is why the 
qq annihilation channel represents a small contribution to the cross sections in the entire 
kinematical domain reached at RHIC and LHC, see e.g. [4]. As a consequence, direct 
photon production is very sensitive to the gluon content of protons and nuclei at these 
colliders [17]. 

At next-to-leading order (NLO), 0{acmCtf), there are altogether 7 types of contributing 
direct subprocesses (see [21]): q{q)g iq{q)g, qiqi iqiql, mi IQjf^iQj^ m 199, 
qi{Qi)qjM9j) 7Qi{Qi)QjiQj)^ Qi{Qi)qi{qi) imiltiqiili), and gg -fqiql. At LHC, qg 
Compton scattering as well as gluon fusion are the processes which dominate the total 
cross section due to the high gluon density in the proton at small values of Bjorken-x. 

In addition to the direct production channels discussed above, prompt photons can 
also be emitted through collinear fragmentation from high-p^, quarks or gluons which are 
produced in primary hard partonic collisions. This is, again schematically, expressible as 

dah,h,^^+x = 5^/^(M2) /j'2(M2) ® daij^k+X'{f^^M^,Mj) ® D^/^iUj). (2.2) 

The perturbatively calculable pieces related to the partonic hard processes, d&ij^k+x', 
are now of the order 0{a1'^"') (n = for LO, and 1 for NLO parton production). The 
collinear divergences associated to this source of photons are resummed and absorbed 
into scale-dependent fragmentation functions (FFs) into photons, D^jj^^z, Mp), where z 
is the fractional momentum over which the last convolution is taken, and Mp is the fac- 
torization scale related to the fragmentation process. The FFs scale asymptotically like 
C(«cm/tts) [22, 23], and thus the convolution of the FFs with the (higher-order) partonic 
cross sections makes the fragmentation contributions of the same order as the direct chan- 
nel. 

We present here results from a complete NLO calculation, which includes both the 
direct and the fragmentation processes at this order. The inclusive calculations performed 
here are based on the INCNLO-package [23, 24]. We use the CTEQ6.6M NLO PDF sets for 
the free proton PDFs [25] and the Bourhis-Fontannaz-Guillet fragmentation functions of 
photons [26]. The renormalization and factorization scales we take to be equal and of the 
order of the photon transverse momentum, /i^ = = M"^ = a x p^,. To investigate the 
uncertainty related to the perturbative nature of the calculation presented, we however vary 
the constant a between 1/2 and 2. Finally, the prompt photon results are also contrasted 
to what is expected for inclusive single-pion production, using the AKK08 set [27] of 
fragmentation functions. 

2.2 Nuclear Parton Distribution Functions 

We estimate the nuclear modifications of the photon spectra at RHIC and LHC using the 

recent NLO releases of nuclear PDFs: nDS [10], EPS09 [13] and HKN07 [11]. All these sets 

p/A 

quantify the scale-dependent ratios between the PDF of a proton inside a nucleus, f- , 
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and that in the unbound proton, /?*, 

These modifications are known to exhibit a rich structure as a function of x. The sup- 
pression at small x < 10~^ is commonly referred to as shadowing while an enhancement 
(anti-shadowing) is predicted around x ~ 10~^. At x > 0.3 the ratio becomes again smaller 
than 1 (EMC-effect), and larger than 1 just below x = 1 because of the Fermi- motion in 
nuclei. 




Figure 1. The nuclear modifications of the gluon PDF in a Pb- nucleus at = 1.69 GeV^ (left) 
and = 100 GcV^ (right) in the EPS09 set (and their uncertainties), nDS set and HKN set. 

The uncertainties related to these modifications are still significant especially in the 
low-Q^ gluon sector, as shown in Fig. 1 (left panel) by the EPS09 uncertainty band on 
the ratio R^^ at = 1.69 GeV^. The small-x uncertainties decrease, however, towards 
higher scales but the large-x uncertainties remain, as is shown by the right panel of Fig. 1. 
In the EPS09 set, the uncertainties of Rf{x,Q^) are encoded through 15 pairs of error 
sets, S^^^ whose gluon modifications are plotted as a light blue band in Fig. 1. The 
propagation of the nPDF uncertainties into a physical quantity X, such as the gluon 
modification in Fig. 1 or prompt photon production considered here, can be obtained by 
squaring the deviations from the central result (which is obtained with the best fit, the set 
5*0) using the following prescription [13]: 

(AX+f « Yl [max{X(5+) - X{So),X{Sr) - X{So),0}]\ 

i 

(AX-f ^Y.[m^x{XiSo) - XiS+),XiSo) - XiS-),0}]\ (2.4) 

i 

For more details about the nPDFs and their uncertainties, see Ref. [13]. 
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3 Probing nPDFs in p—A collisions 

3.1 Feasibility of direct photon measurements 

Before discussing the nPDF corrections to the prompt photon yield in nuclear collisions, 
we first briefly discuss the expected rates and hence the statistical accuracy at RHIC and 
LiiC. In Fig. 2 we plot the absolute NLO prompt photon p^-spectra at mid-rapidity in 
(minimum-bias) d-Au collisions at RHIC and p-Ph collisions at the LHC, together with 
their scale uncertainty (band). Assuming that it would be necessary to reach a 6 = 10% 
statistical accuracy in order to measure meaningfully the nuclear effects in prompt photon 
production, one would need about M = w 100 events in a {p^,, rapidity)-bin. 




10 20 30 40 50 60 70 80 90 100 



Pt (GeV) 



Figure 2. da/dpj,dy inclusive photon production -spectra in minimum-bias d-Au collisions at 
RHIC (ys;;;; = 200 CeV) and p-Pb collisions at LHC = 8.8 TeV). The band (not visible at 

LHC energy) corresponds to a scale variation from n = M = Mp = [p^/^; 2^^]. 

The foreseen luminosity in d-Au collisions at RHIC-II is C^^'^'^ = 0.74 pb~^ assuming 
12 weeks of ion runs per year [28]. Requiring M = 100 events/GeV, the minimal cross 
section reads 

which would be reached at Pj. ~ 35 GeV, see Fig. 2. At LHC such a statistical accuracy 
would correspond to a cross section 

_ = _ pb/GeV . 10^' pb/GeV, (3.2) 

assuming a rather conservative C = 10^^ cm~^s~^ leading to a yearly integrated integrated 
luminosity^ Z:'"^* = 0.1 pb~^ [29]. From Fig. 2, this precision is achieved for transverse 



^We assume here that the LHC will run one month per year in the ion mode, which is taken to be 
At = 10^ s by convention. 
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momenta up to p^, ~ 100 GeV in one year. These considerations set the Pj,-windows for 
the LHC and RHIC in which we perform our computations. Better precision could only be 
achieved of course in a more hmited Pj. range. A more satisfactory 5 = 3% precision would 
be reached for transverse momenta less than Pj, ~ 25 GeV at RHIC and Pj, ~ 60 GeV at 
the LHC. 

3.2 RHIC 

Let us now move to the prediction of the nuclear modification ratio, or "quenching factor" , 
in d-Au collisions at RHIC, defined as 

„7 ^ da/dp^ ((i + Au 7 + X) 
'^A'^ ^ 2 X 197 X da/dpy (p + p^ 7 + X)' 

Note that the quenching factor (3.3) is not normalized to one when no nuclear modifica- 
tions in the parton densities are assumed. The reason comes from an "isospin" effect as 
the density of up quarks - to which photons mostly couple - in the deuterium and the 
gold nucleus is smaller than that in a proton because of the presence of neutrons inside 
nuclei [30]. These corrections should be most pronounced whenever the valence quark sec- 
tor of the nuclei is probed, that is at large xi/x2 and thus large Pj,. On the contrary, no 
significant isospin effect is expected in charge-averaged or neutral hadron production which 
do not involve QED couplings. 




5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 
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Figure 3. Nuclear modification ratio RJj^^ of inclusive photon production at \y\ < 0.35 in d-Au 
collisions at y/s^ = 200 GeV using EPS09 nPDFs, in comparison to (i) uDS and HKN sets (left), 
and (ii) the pion case, R^j^^ (right). 

The predictions for -R^^^ at mid-rapidity and ^/s^ = 200 GeV are shown in Fig. 3 
(left) in the transverse momentum range Pj. = 5-50 GeV. The central EPS09 prediction is 
shown as a dashed black line, the light blue band corresponding to its uncertainty range, 
while the nDS (resp. HKN) prediction is denoted by the dotted brown (resp. dash-dotted 
green) curve. For separating the genuine nuclear effects from the sheer isospin effects, the 
solid red line indicates the calculation with no nuclear modifications in the PDFs (i.e. just 
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the free proton PDFs are used) - it also serves as a comparison baseline for the nuclear 
effects obtained using the various nPDF sets. 

The relative difference between EPS09 with respect to the baseline (red curve) predic- 
tion follows roughly the shape of the nuclear modifications in Fig. 1. Below Pj. < 10 GeV 
photon production is sensitive to anti-shad owing corrections {x2 ~ 2pj,/^s^^ ~ 0.1) while 
it is suppressed at larger momenta due to the EMC effect. This is also the case for pion 
production for which we only show the EPS09 prediction on as a blue band in the 

right panel of Fig. 3, except that in this case the suppression is practically free of isospin 
corrections, as mentioned above. 
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Figure 4. Same as Fig. 3 at y = 3. 



As noted in section 2.1, prompt photon production at forward rapidities, y > 0, is 
more sensitive to lower momentum fractions in a target nucleus, X2 oc exp(— y), than the 
mid-rapidity production, and therefore it probes efficiently the gluon nPDFs at small x. 
We demonstrate this in Fig. 4 (left) by plotting the analogue of Fig. 3 (left) at rapidity 
y = 3 and in a narrower Pj,-range due to phase-space restriction. As expected, -RJ^^ at 
small Pj, falls now below the baseline (red) curve, signaling that this kinematical domain 
lies in the small-x shadowing region [p^ = 5 GeV and y = 3 corresponds to a typical 
value X2 ~ 2p^e~y / ^s^j^ ~ 2 x 10^'^). From these figures, however, one may see that the 
predicted differences due to nPDF effects are not very large - of the order of 10% or less. 
For completeness, the quenching factor of single-pion production, R^^^, at y = 3 is also 
plotted together with RJj^^ in Fig. 4 (right). The behaviour of is again roughly similar 
to RJj^^ except for the isospin corrections. 

It might be surprising at first glance to observe a strong isospin effect in the prompt 
photon channel at forward rapidity since quark distributions in a proton and in a neutron 
are symmetric at small values of x: vP{x) = (i"(x) ~ dP{x) = u^{x). As a matter of fact, 
the strong isospin corrections visible in Fig. 4 actually come from the deuteron projectile 
which is probed at large xi oc exp{+y). As we shall see in the next section, unlike in 
the RHIC case, the isospin effect vanishes at forward rapidity at the LHC since the same 
projectile is used in p-p and p-Fh collisions. 
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Figure 5. Left: Prompt photon quenching factor RJj^^ as a function of in different rapidities: 
y ^ (sohd), y = 1.6 (dotted) and y = 2.4 (dash-dotted) with EPS09 nPDF corrections. Right: 

(proton PDFs). 
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The RHIC photon data wih be useful to learn about the behaviour of the gluon nPDFs 
from small to large x. In order to show the different kinematical regions probed at forward 
rapidities, -RJ^^ is plotted in Fig. 5 (left) as a function of p^, for three rapidity bins {y = 0, 
y = 1.6, and y = 2.4) which should be accessible by the PHENIX experiment at RHIC. 
The predictions are obtained assuming the EPS09 nPDFs. The EPS09 -R^^^ ratio looks 
remarkably flat at all rapidity bins, even though the normalization differs: the larger 
the rapidity, the stronger the shadowing and the smaller the quenching factor. This flat 
behaviour is however somewhat accidental as it results from the interplay of isospin and 
nPDF effects. As can be seen in Fig. 5 (right) where the quenching factor is normalized 
with what is expected assuming isospin effects only, the behaviour is very different^. At 
mid-rapidity, prompt photon production is sensitive to anti-shadowing effects in this p^,- 
range (see also Fig. 3) whereas the crossing region between anti-shadowing and shadowing 
is probed at y = 1.6 and y = 2.4 around Pj, = 5-6 GeV. 

3.3 LHC 

In this section we give predictions for p-Pb collisions at the LHC. The nuclear modification 
factor. 



da/dpj. (p -h Pb 7 X) 



2QSxda/dp^{p + p^-i + Xy ^'^'^^ 
is plotted for prompt photon production at mid-rapidity in Fig. 6. Because of the larger 



center-of-mass energy. 



8800 GeV, the typical values of X2 oc 1/ probed in the 



nuclear target are much smaller than at RHIC. Remarkably, the differences in the results 
obtained with the different nPDF sets, are now more pronounced. 



^This would correspond roughly to the production ratio oc d?a [d + Au 
in which isospin corrections are very small. 



■7 + X) /dV(d + d^7 + X) 
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In the EPS09 set (light blue band), prompt photon production is enhanced with respect 
to the isospin effects (solid red line) above p^, = 20-40 GeV due to the anti-shadowing. On 
the contrary, the transition from the shadowing to the anti-shadowing regions occurs only 
above Pj, ~ 120 GeV when using nDS (brown short-dashed) and HKN (green dot-dashed) 
nPDF sets. Like at RHIC, it is particularly interesting to note that the predictions using the 
latter two sets fall outside the EPS09 uncertainty band. As a consequence, this observable 
- which should be easily accessible at LHC - will allow one to set tighter constraints to 
the nuclear gluon densities. The EMC-effect (appearing at X2 > 0.1) which can be probed 
easily at RHIC would become visible at the LHC only at very large p^, Pj. > 500 GeV at 
y = 0, or at negative rapidities, say, at y < —2.2 for p^ = 50 GeV. 
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Figure 6. Nuclear modification ratio RJj^^ of inclusive photon production at \y\ < 0.5 in p-Pb 
collisions at y/s^^ = 8.8 TeV using the EPS09 nPDF set, in comparison (left) to nDS and HKN 
nPDFs and (right) to the pion production channel. 

The smallest X2 values can be reached at LHC (because of the high c.m.s. energy), at 
forward rapidity and at not too large values of p^, . This kinematic region is therefore ideal 
for probing efficiently the shadowing region. This can be seen in Fig. 7 (left) where the 
photon suppression at y = 3 is plotted as a function of p^,. As discussed in the previous 
section, no isospin corrections (solid red line) are expected since at small values of X2 the 
quark distributions in protons and neutrons are identical. Without a surprise, the photon 
quenching factor is below 1 because of the shadowing in all the nPDF sets studied here. 
The large uncertainty band of the EPS09 predictions in Fig. 7 reflects that of R^^ shown 
in Fig. 1. Therefore, high-precision measurements at p^, ~ 10 GeV and at y = 3 in p- 
Pb collisions at LHC would bring significant further constraints for the nPDF global fit 
analyses. 

Going to the backward direction at LHC would allow one to probe nPDFs in an x-region 
similar to that reached with mid-rapidity photon production at RHIC energy. As discussed 
in [17], measuring prompt photon production in p-A collisions at negative rapidity would 
allow one to access the quark nPDFs at large values of X2, R"^'^°{x^) ~ i?*^ {x^ exp(— y)), 
while mid-rapidity photon production is sensitive to both the quark and gluon nPDFs, 
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■^pA (^-l) ~ ((-^Fj^^i-^ + Fig- 7 (right) shows the expected photon suppression 

in the rapidity bin y = —3 in the range Pj. = 10-100 GeV. This range would correspond at 
RHIC, at similar values of X2, to transverse momenta Pj, = O (-^/Snn^™^/ V'^nn^^*^ ^ ~ 
5-50 GeV at mid-rapidity. As can be seen from the comparison of Fig. 7 (right) with Fig. 3 
(left), the expected photon suppression is rather similar. Remarkably, the spread of the 
EPS09 theoretical predictions proves narrower at the LHC than at RHIC, reflecting the 
fact that quark nPDFs are much better constrained than the gluon nPDFs at large values 
of X [13]3. 
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Figure 7. Nuclear production ratio of inclusive photon production at ?; = 3 (left) and y = —'i 
(right) in p-Pb collisions at ^/s^ ~ 8.8 TeV using the EPS09 nPDF set, in comparison to the uDS 
and HKN nPDF sets. 



3.4 Scale dependence 

As discussed in section 2.1, collinearly factorized pQCD cross sections depend on the renor- 
malization and factorization scales, which all are of the order of the photon p^. In the 
absolute cross sections the sensitivity to these scales reflects the uncertainty which results 
from terminating the perturbation series at a certain order, and thus neglecting the higher- 
order corrections (here, NNLO and beyond). In the nuclear modification ratio RJ^^, such 
scale uncertainties should nevertheless largely just cancel out. However, since the nPDF 
corrections Rf{x,M^) (Eq. (2.3)) do depend on the factorization scale M, also RJj^^ may 
exhibit some dependence on M. To quantify this theoretical uncertainty, the nuclear mod- 
ification ratio of prompt photon production at mid-rapidity has been computed using the 
EPS09 set and varying all scales, /u = M = Mp from Py/2 to 2pj,, as was done also in the 
calculation of the absolute cross section in section 3.1. As shown in Fig. 8 at the RHIC 
and LHC energies, the predictions show very little scale dependence. More importantly, 
the scale dependence proves much smaller than the current uncertainties in the nuclear 
modifications of the PDFs; see the band in Fig. 2. 

^In addition, larger scales are probed at the LHC, Q'^Iljjq ~ 4 x Q'^Ij^jjjq- This is however a rather 
moderate effect since the EPS09 gluon nPDF ratios do not exhibit a strong Q'^-dependence at large x [13]. 
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Figure 8. Scale dependence of the nuclear modification ratio for mid-rapidity prompt photon 
production at RHIC (lower curve) and LHC (upper curve). Calculations are performed using the 
central EPS09 nPDF set and varying all scales by a factor of two with respect to the photon 
transverse momentum; see the text for details. 

4 Baseline pQCD predictions in heavy-ion collisions 
4.1 Inclusive prompt photon production 

As can be seen in Eqs. (2.1) and (2.2), nuclear corrections are expected to be much more 
pronounced in nucleus-nucleus collisions than in proton-nucleus collisions, very roughly 
Raa ~ RpA ~ i^f)'^ ™ minimum-bias collisions. Therefore, it is particularly useful to 
investigate how prompt photon production in heavy-ion collisions would show the presence 
of a strongly interacting QCD-medium and possibly also constrain nuclear parton densities. 
This would also allow for more reliable estimates of other hard processes like jet and large- 
Pj, hadron production, which, on the basis of RHIC data (see e.g. [31] for a review) and 
the first LHC measurements [32, 33], are generally believed to be affected by energy losses 
of hard partons in the formed quark-gluon plasma (QGP). Another motivation to make 
the baseline pQCD predictions for the A-A collisions here is that the Pb-Pb programme 
at the LHC has already started. 

The quenching factors in Au-Au collisions at RHIC (left) and Pb-Pb collisions at the 
LHC (right) are shown in Fig. 9. As expected, the shapes are similar but the nuclear 
effects are stronger than in the respective p-A collisions, as the comparison with Fig. 3 
(left) and Fig. 6 (left) shows. Note moreover that the nPDF corrections might be even 
more pronounced in central collisions as compared to these minimum-bias predictions, 
which are averaged over all centralities. A study of the centrality dependence of -R^^^^ 
would require a modeling of the spatial dependences of the nPDFs, which is outside the 
scope of the present study. The PHENIX collaboration at RHIC has reported on high-pj, 
measurements of -R^uAu central collisions [34] (see Fig. 2 there). However, since these 
data are still preliminary and since our present setup is consistent only with minimum-bias 



- 11 - 



collisions, we do not make a detailed comparison with these data in this paper. 



1.6 
1.4 

< 1.2 



1 — f — 1 — f — , — f — , — I — , — I — , — f — , — f — 1 — f — , — f — 

lykO.35 proton PDF- 

1/2 ^ , , EPS09 

s =200 GeV --. nDS 

HKN 



0^ 



10 15 



20 25 30 35 

Pt [GeV] 



1.2 
1.1 
1.0 
0.9 
0.8 
0.7 



40 45 50 



1 1 1 1 1 1 1 1 
lykO.5 

1/7 

■ s =5.5TeV 


1 1 1 1 1 1 1 1 1 
proton PDF 

EPS09 

nDS 


— 


HKN 














20 40 60 



80 100 120 140 160 180 200 

Pt [GeV] 



Figure 9. Left: Nuclear production ratio i?2uAu '^'^ inclusive photon production in minimum-bias 
Au-Au collisions at =■ 200 GeV (|y| < 0.35) using EPS09, nDS and HKN nuclear parton 

densities. Right: same in Pb-Pb collisions at Js^^ = 5.5 TeV (|?;| < 0.5). 



The plots in Fig. 9 thus set the stage for the conclusions to be made from the forthcom- 
ing data: agreement with the pQCD baseline in p-A collisions but not in A-A collisions 
would be a clear hint that photons are actually sensitive to hot medium effects. 

4.2 Quenching of the fragmentation component 

It is not clear yet to what extent the production of prompt photons is affected by the 
presence of the dense medium in heavy-ion collisions. On the one hand, it has been sug- 
gested that the parton multiple scattering might enhance the photon yield either through 
jet-to-photon conversion [35] or medium-induced photon bremsstrahlung [36]. On the other 
hand, photon production might be suppressed due to the quenching of the fragmentation 
component as it happens for inclusive hadron production [30]. 

At large enough p^,, say at p^, > 5 GeV at RHIC and possibly slightly above at LHC, 
one could nevertheless expect the suppression due to the quenching of the fragmentation 
component to be larger than the possible enhancement (s) caused by the medium. There- 
fore, we present in Fig. 10 a set of curves where we have hy hand simply downscaled the 
fragmentation component by a factor of five (assuming in the calculation the central EPS09 
set, for the illustration) which very roughly corresponds to the observed maximum suppres- 
sion of large-pj, hadrons at RHIC and LHC. These curves should therefore serve as a lower 
limit for prompt photon production in A-A collisions. It should, however, be noted that 
these additional curves are, strictly speaking, not physical quantities since the separation 
between direct and fragmentation components is ambiguous beyond the Born level and 
scale-dependent. As before, the calculations are performed for different scale choices from 
Pj,/2 to 2pj, corresponding to the band in Fig. 10. We also reproduce in this figure the 
quenching factor obtained with the EPS09 nPDF sets (with its EPS09 uncertainty band) 
as shown in Fig. 9, i.e. without any rescaling of fragmentation component, for comparison. 
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Figure 10. Upper bands: Same as Fig. 9 using the EPS09 nPDF sets. Lower bands: Nuclear 
production ratio R\j^ using the central EPS09 nPDF set and after downscaling the fragmentation 
component by a factor of five. The band corresponds to a scale variation from Pj,/2 to . 



The prompt photon suppression is most pronounced at small p^., where the fragmen- 
tation component is the largest. As Pj, gets larger, the direct photon contribution becomes 
more important, and therefore the expected photon quenching weakens. At RHIC, the 
photon quenching becomes comparable to the prediction with no energy loss effects above 
Pj, > 25 GeV. The RHIC predictions are qualitatively similar to the calculation performed 
in [30]. At LHC, the direct photon contribution to the total cross section is somewhat less 
than at RHIC. As a consequence, energy loss effects are more pronounced and might be 
visible in the entire p^-range. 



5 Summary 

Prompt photon production in p-A and d-A collisions is among the best observables for 
probing parton distributions in nuclei, especially the nuclear gluon distributions, which are 
presently still poorly known. In this paper, we have computed the nuclear modification 
factors of single inclusive prompt photon production at NLO accuracy using the recent 
nPDF sets EPS09, nDS, and HKN. Calculations were carried out for d-Au collisions at 
RHIC and p-Ph collisions at the LHC, for their nominal energies. The results are conve- 
niently complementary: At mid-rapidities, prompt photon production at RHIC is mostly 
sensitive to the gluon anti-shad owing and the EMC effect, while at the LHC these photons 
probe the small-x shadowing region. Complementary information on the gluon shadowing 
region can be obtained also from the forward-rapidity prompt photons both at RHIC and 
the LHC. Remarkably, by combining the future data at both colliders, one should be able 
to discriminate among the different existing sets of nPDFs. 

On top of constraining the nuclear PDFs, inclusive prompt photon production can 
also be used as a "calibration measurement" in A-A collisions where QGP formation is 
expected to quench the production of large-p^ partons, and thereby also the fragmentation 
component of prompt photon production. We have reported two pQCD baselines for 
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this purpose: one by including the nuclear effects to the PDFs but no medium-effects to 
prompt photon production, and one by using the nPDFs and suppressing the fragmentation 
component by a factor of five. 

The photons produced by medium-induced bremsstrahlung or through fragmentation 
processes are likely to be produced alongside a significant hadronic activity in their vicinity. 
This is also true - perhaps to a lesser extent - with the jet-photon conversion process, since 
the initial parton will radiate before "converting" into a photon though a rescattering in 
the QGP. Therefore, a way to reduce their contributions - and therefore medium effects - 
would be to trigger on isolated photons, i.e. using isolation criteria around each photon 
candidate. In heavy-ion collisions, the high-multiplicity would make this measurement ex- 
tremely challenging, hopefully made possible using recent advances developed towards the 
reconstruction of jets in heavy-ion collisions [37]. Predictions of isolated photon produc- 
tion in nuclear collisions go beyond the scope of the present paper and will be presented 
elsewhere. 
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